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Abstract. Observational advances over the last decade have left little doubt that neu- 
tron stars received a large kick velocity (of order a few hundred to a thousand kms ) 
at birth. The physical origin of the kicks and the related supernova asymmetry is one of 
the central unsolved mysteries of supernova research. We review the physics of different 
kick mechanisms, including hydrodynamically driven, neutrino - magnetic field driven, 
and electromagnetically driven kicks. The viabilities of the different kick mechanisms 
are directly related to the other key parameters characterizing nascent neutron stars, 
such as the initial magnetic field and the initial spin. Recent observational constraints 
on kick mechanisms are also discussed. 

1 Evidence for Neutron Star Kicks and Supernova 
Asymmetry 

It has long been recognized that neutron stars (NSs) have space velocities much 
greater (by about an order of magnitude) than their progenitors', (e.g., Gunn & 
Ostriker 1970). A natural explanation for such high velocities is that supernova 
explosions are asymmetric, and provide kicks to the nascent NSs. In the last 
few years, evidence for NS kicks and supernova asymmetry has become much 
stronger. The observational facts and considerations that support (or even re- 
quire) NS kicks fall into three categories: 

(1) Large NS Velocities (^> the progenitors' velocities ~ 30 kms -1 ): 

• Recent studies of pulsar proper motion give 200—500 kins" 1 as the mean 3D 
velocity of NSs at birth (e.g., Lyne and Lorimer 1994; Lorimer et al. 1997; Hansen 
& Phinney 1997; Cordes & Chernoff 1998), with possibly a significant population 
having velocities greater than 1000 kms -1 . While velocity of ~ 100 kms -1 may 
in principle come from binary breakup in a supernova (without kick), higher 
velocities would require exceedingly tight presupernova binary. Statistical anal- 
ysis seems to favor a bimodal pulsar velocity distribution, with peaks around 
100 kms -1 and 500 kms -1 (Arzoumanian et al. 2001; see also Hansen & Phin- 
ney 1997; Cordes & Chernoff 1998). 

• Direct evidence for pulsar velocities >1000 km s _1 has come from obser- 
vations of the bow shock produced by the Guitar Nebula pulsar (B2224+65) in 
the interstellar medium (Cordes, Romani & Lundgren 1993). 

• The studies of neutron star - supernova remnant associations have, in many 
cases, indicated large NS velocities (e.g., Frail et al. 1994), although identifying 
the association can be tricky sometimes (e.g. Kaspi 1999; Gaensler 2000). Of spe- 
cial interest is the recent studies of magnetar-SNR associations: the SGR 0526-66 
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- N49 association, implying V± ~ 2900 (3 kyr/t) km s 1 , and the possible associ- 
ation of SGR 1900+14 with G42. 8+0.6, implying V± ~ 1800 (10 kyr/i) kms" 1 . 
(However, the proper motion of SGR 1806-20 may be as small as 100 km s _1 , 
and AXP 1E2259+586, AXJ 1845-0258, and AXP 1E1841-045 lie close to the 
centers of their respective remnants, CTB 109, G29. 6+0.1, and Kes 73) (see 
Gaensler 2000). 

(2) Characteristics of NS Binaries (Individual Systems and Popula- 
tions): While large space velocities can in principle be accounted for by binary 
break-up (as originally suggested by Gott et al. 1970; see Iben & Tutukov 1996), 
many observed characteristics of NS binaries demonstrate that binary break-up 
can not be solely responsible for pulsar velocities, and that kicks are required 
(see also Tauris & van den Heuvel 2000). Examples include: 

• The detection of geodetic precession in binary pulsar PSR 1913+16 implies 
that the pulsar's spin is misaligned with the orbital angular momentum; this can 
result from the aligned pulsar-He star progenitor only if the explosion of the He 
star gave a kick to the NS that misalign the orbit (Cordes et al. 1990; Kramer 
1998; Wex et al. 1999). 

• The spin-orbit misalignment in PSR J0045-7319/B-star binary, as mani- 
fested by the orbital plane precession (Kaspi et al. 1996; Lai et al. 1995) and fast 
orbital decay (which indicates retrograde rotation of the B star with respect to 
the orbit; Lai 1996a; Kumar & Quataert 1997) require that the NS received a 
kick at birth (see Lai 1996b). 

• The observed system radial velocity (430 km s -1 ) of X-ray binary Circinus 
X-l requires Vidck > 500 km s" 1 (Tauris et al. 1999). 

• High eccentricities of Be/X-ray binaries cannot be explained without kicks 
(Verbunt & van den Heuvel 1995). 

• Evolutionary studies of NS binary population (in particular the double NS 
systems) imply the existence of pulsar kicks (e.g., Deway & Cordes 1987; Fryer 
& Kalogera 1997; Fryer et al. 1998). 

(3) Observations of SNe and SNRs: 

There are many direct observations of nearby supernovae (e.g., spectropo- 
larimetry: Wang et al. 2000, Leonard et al. 2000; X-ray and gamma-ray observa- 
tions and emission line profiles of SN1987A: McCray 1993, Utrobin et al. 1995) 
and supernova remnants (e.g., Morse, Winkler & Kirshner 1995; Aschenbach et 
al. 1995) which support the notion that supernova explosions are not spherically 
symmetric. 

Finally it is of interest to note that recent study of the past association of 
the runaway star C Oph with PSR J1932+1059 (Hoogerwerf et al. 2000) or with 
RX 185635-3754 (Walter 2000) also implies a kick to the NS. 

2 The Problem of Core-Collapse Supernovae and 
Neutron Star Kicks 

The current paradigm for core-collapse supernovae leading to NS formation is 
that these supernovae are neutrino-driven (see Bethe 1990, Burrows 2000, Janka 
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2000 for recent review): As the central core of a massive star collapses to nuclear 
density, it rebounds and sends off a shock wave, leaving behind a proto-neutron 
star. The shock stalls at several 100's km because of neutrino loss and nuclear 
dissociation in the shock. A fraction of the neutrinos emitted from the proto- 
neutron star get absorbed by nucleons behind the shock, thus reviving the shock, 
leading to an explosion on the timescale several 100's ms — This is the so-called 
"Delayed Mechanism". However, ID simulations with detailed neutrino trans- 
port seem to indicate that neutrino heating of the stalled shock, by itself, does not 
lead to an explosion or produce the observed supernova energetics (see Rampp 
& Janka 2000) . It has been argued that neutrino-driven convection in the proto- 
neutron star (which tends to increase the neutrino flux) and that in the shocked 
mantle (which tends to increase the neutrino heating efficiency) are central to the 
explosion mechanism, although there is no consensus on the robustness of these 
convections (e.g., Bethe 1990; Herant et al. 1994; Burrows et al. 1995; Janka & 
Muller 1996; Mezzacappa et al. 1998). What is even more uncertain is the role 
of rotation and magnetic field on the explosion (see Monchmeyer et al. 1991; 
Rampp, Muller & Ruffert 1998; Khokhlov et al. 1999; Fryer & Heger 2000 for 
simulations of collapse/explosion with rotation, and Thompson & Duncan 1993 
and Thompson 2000a for discussion of possible dynamo processes and magnetic 
effects). 

It is clear that despite decades of theoretical investigations, our understand- 
ing of the physical mechanisms of core-collapse supernovae remains significantly 
incomplete. The prevalence of neutron star kicks poses a significant mystery, and 
indicates that large-scale, global deviation from spherical symmetry is an impor- 
tant ingredient in our understanding of core-collapse supernovae (see Burrows 
2000). 

In the following sections, we review different classes of physical mechanisms 
for generating NS kicks (§§3-5), and then discuss possible observational con- 
straints and astrophysical implications (§6). 

3 Hydrodynamically Driven Kicks 

The collapsed stellar core and its surrounding mantle are susceptible to a variety 
of hydrodynamical (convective) instabilities (e.g., Herant et al. 1994; Burrows 
et al. 1995; Janka & Muller 1996; Keil et al. 1996; Mezzacappa et al. 1998). 
It is natural to expect that the asymmetries in the density, temperature and 
velocity distributions associated with the instabilities can lead to asymmetric 
matter ejection and/or asymmetric neutrino emission. Numerical simulations, 
however, indicate that the local, post-collapse instabilities are not adequate to 
account for kick velocities >100 km s _1 (Janka & Muller 1994; Burrows & Hayes 
1996; Janka 1998; Keil 1998) — These simulations were done in 2D, and it is 
expected that the flow will be smoother on large scale in 3D simulations, and 
the resulting kick velocity will be even smaller. 

There is now a consensus that global asymmetric perturbations in presuper- 
nova cores are required to produce the observed kicks hydrodynamically (Gol- 
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dreich, Lai & Sahrling 1996; Burrows & Hayes 1996). Numerical simulations 
by Burrows & Hayes (1996) demonstrate that if the precollapse core is mildly 
asymmetric, the newly formed NS can receive a kick velocity comparable to 
the observed values. (In one simulation, the density of the collapsing core ex- 
terior to 0.9Mq and within 20° of the pole is artificially reduced by 20%, and 
the resulting kick is about 500 km s _1 .) Asymmetric motion of the exploding 
material (since the shock tends to propagate more "easily" through the low- 
density region) dominates the kick, although there is also contribution (about 
10 — 20%) from asymmetric neutrino emission. The magnitude of kick velocity 
is proportional to the degree of initial asymmetry in the imploding core. Thus 
the important question is: What is the origin of the initial asymmetry? 

3.1 Presupernova Perturbations 

Goldreich et al. (1996) suggested that overstable g-mode oscillations in the pre- 
supernova core may provide a natural seed for the initial asymmetry. These 
overstable g-modes arise as follows. A few hours prior to core collapse, a massive 
star (M > 8Mq) has gone through a successive stages of nuclear burning, and 
attained a configuration with a degenerate iron core overlaid by an "onion skin" 
mantle of lighter elements. The rapidly growing iron core is encased in and fed 
by shells of burning silicon and oxygen, and the entire assemblage is surrounded 
by a thick convection zone. The nearly isothermal core is stably stratified and 
supports internal gravity waves. These waves cannot propagate in the unstably 
stratified convection zone, hence they are trapped and give rise to core g-modes 
in which the core oscillates with respect to the outer parts of the star. The over- 
stability of the g-mode is due to the "e-mechanism" with driving provided by 
temperature sensitive nuclear burning in Si and O shells surrounding the core 
before it implodes. It is simplest to see this by considering a I = 1 mode: If we 
perturb the core to the right, the right-hand-side of the shell will be compressed, 
resulting in an increase in temperature; Since the shell nuclear burning rate 
depends sensitively on temperature (power-law index <~ 47 for Si burning and 
~ 33 for O burning), the nuclear burning is greatly enhanced; this generates a 
large local pressure, pushing the core back to the left. The result is an oscillating 
g-mode with increasing amplitude. 

The main damping mechanism comes from the leakage of mode energy. The 
local (WKB) dispersion relation for nonradial waves is 

kl = {Jcl)-\J-tf){J-N% (1) 

where k r is the radial wavenumber, L\ — + l)c s /r (c s is the sound speed) 
and N are the acoustic cut-off (Lamb) frequency and the Brunt- Vaisala, fre- 
quency, respectively. Since acoustic waves whose frequencies lie above the acous- 
tic cutoff can propagate through convective regions, each core g-mode will couple 
to an outgoing acoustic wave, which drains energy from the core g-modes (see 
Fig. [j]). This leakage of mode energy can be handled with an outgoing propaga- 
tion boundary condition in the mode calculation. Also, neutrino cooling tends to 
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Fig. 1. Propagation diagram computed for a 15Mq presupernova model of Weaver and 
Woosley (1993). The solid curve shows TV 2 , where N is the Brunt- Vaisala frequency; 
the dashed curves show Lf, where L; is the acoustic cutoff frequency, with 1 = 1, 2, 3. 
The spikes in iV 2 result from discontinuities in entropy and composition. The iron core 
boundary is located at 1.3M , the mass-cut at 1.42Mq. Convective regions correspond 
to N = 0. Gravity modes (with mode frequency ui) propagate in regions where ui < N 
and lu < Li, while pressure modes propagate in regions where w > N and w > Li . Note 
that a g-mode trapped in the core can lose energy by penetrating the evanescent zones 
and turning into an outgoing acoustic wave (see the horizontal line). Also note that 
g-modes with higher n (the radial order) and I (the angular degree) are better trapped 
in the core than those with lower n and I. 



damp the mode. Since the nuclear energy generation rate depends more sensi- 
tively on temperature than pair neutrino emission (power law index ~ 9) , cooling 
is never comparable to nuclear heating locally. Instead, thermal balance is me- 
diated by the convective transport of energy from the shells, where the rate of 
nuclear energy generation exceeds that of neutrino energy emission, to the cooler 
surroundings where the bulk of the neutrino emission takes place. Calculations 
(based on the 15M and 25M presupernova models of Weaver & Woosley 1993) 
indicate that a large number of g-modes are overstablc, although for low-order 
modes (small I and n) the results depend sensitively on the detailed structure 
and burning rates of the presupernova models. The typical mode periods are 
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>1 s, the growth time ^10 — 50 s, and the lifetime of the Si shell burning is ~ 
hours (Lai & Goldreich 2000b, in preparation). 

Our tentative conclusion is that overstable g-modes can potentially grow to 
large amplitudes prior to core implosion, although a complete understanding of 
the global pre-collapse asymmetries is probably out of reach at present, given the 
various uncertainties in the presupernova models. For example, the O-Si burning 
shell is highly convective, with convective speed reaching 1/4 of the sound speed, 
and hydrodynamical simulation may be needed to properly modeled such con- 
vection zones (see Bazan & Arnett 1998, Asida & Arnett 2000). Alternatively, it 
has been suggested that the convection itself may provide the seed of asymmetry 
in the presupernova core (Bazan & Arnett 1998), although it is not clear whether 
the perturbations have sufficiently large scales to be relevant to supernova kicks. 

3.2 Amplification of Perturbation During Core Collapse 

Core collapse proceeds in a self-similar fashion, with the inner core shrinking 
subsonically and the outer core falling supcrsonically at about half free-fall speed 
(Goldreich & Weber 1980; Yahil 1983). The inner core is stable to non-radial 
perturbations because of the significant role played by pressure in its subsonic 
collapse (Goldreich & Weber 1980). Pressure is less important in the outer region, 
making it more susceptible to large scale instability. A recent stability analysis 
of Yahil's self-similar collapse solution (which is based on Newtonian theory 
and a polytropic equation of state P oc p r , with r ~ 1.3) does not reveal any 
unstable global mode before the proto-neutron star forms (Hanawa & Matsumoto 
2000; Lai 2000). However, during the subsequent accretion of the outer core 
(involving 15% of the core mass) and envelope onto the proto-neutron star, 
nonspherical Lagrangian perturbations can grow according to Ap/p oc r -1 / 2 
(independent of I) or even Ap/p oc r^ 1 (for I = 1 when the central collapsed 
object is displaced from the origin of the converging flow) (Lai & Goldreich 
2000) The asymmetric density perturbations seeded in the presupernova star, 
especially those in the outer region of the iron core, are therefore amplified (by a 
factor of 5-10) during collapse. The enhanced asymmetric density perturbation 
may lead to asymmetric shock propagation and breakout, which then give rise 
to asymmetry in the explosion and a kick velocity to the NS (see Burrows & 
Hayes 1996). 

4 Neutrino — Magnetic Field Driven Kicks 

The second class of kick mechanisms rely on asymmetric neutrino emission in- 
duced by strong magnetic fields. Since 99% of the NS binding energy (a few times 
10 53 erg) is released in neutrinos, tapping the neutrino energy would appear to 
be an efficient means to kick the newly-formed NS. The fractional asymmetry 
a in the radiated neutrino energy required to generate a kick velocity Vkick is 
a = MV kick c/E tot (= 0.028 for V kick = 1000 km s _1 , NS mass M = 1.4 M Q and 
total neutrino energy radiated E tot = 3 x 10 53 erg). 
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4.1 Effect of Parity Violation 

Because weak interaction is parity violating, the neutrino opacities and emissiv- 
ities in a magnetized nuclear medium depend asymmetrically on the directions 
of neutrino momenta with respect to the magnetic field, and this can give rise to 
asymmetric neutrino emission from the proto-neutron star. Chugai (1984) (who 
gave an incorrect expression for the electron polarization in the relativistic, de- 
generate regime) and Vilenkin (1995) considered neutrino-electron scattering, 
but this is less important than neutrino-nucleon scattering in determining neu- 
trino transport in proto-neutron stars. Dorofeev et al. (1985) considered neutrino 
emission by Urea processes, but failed to recognize that in the bulk interior of 
the star the asymmetry in neutrino emission is cancelled by that associated with 
neutrino absorption (Lai & Qian 1998a). 

Horowitz & Li (1998) suggested that large asymmetries in the neutrino flux 
could result from the cumulative effect of multiple scatterings of neutrinos by 
slightly polarized nucleons (see also Lai & Qian 1998a; Janka 1998). However, 
it can be shown that, although the scattering cross-section is asymmetric with 
respect to the magnetic field for individual neutrinos, detailed balance requires 
that there be no cumulative effect associated with multiple scatterings in the 
bulk interior of the star where thermal equilibrium is maintained to a good 
approximation (Arras & Lai 1999a; see also Kusenko et al. 1998). For a given 
neutrino species, there is a drift flux of neutrinos along the magnetic field in 
addition to the usual diffusive flux. This drift flux depends on the deviation of 
the neutrino distribution function from thermal equilibrium. Thus asymmetric 
neutrino flux can be generated in the outer region of the proto-neutron star 
(i.e., above the neutrino-matter decoupling layer, but below the neutrinosphere) 
where the neutrino distribution deviates significantly from thermal equilibrium. 
While the drift flux associated with v^s and i^v's is exactly canceled by that 
associated with D^s and P T 's, there is a net drift flux due to v^s and D e 's. 
Arras & Lai (1999b) found that the asymmetry parameter for the v e -v e flux is 
dominated for low energy neutrinos (<15 MeV) by the effect of ground (Landau) 
state electrons in the absorption opacity, e a b s — 0.6Bi 5 (E u /l MeV) -2 , where 
Bi 5 = B/(10 15 G), and for high energy neutrinos by nucleon polarization (e <~ 
fi m B/T). Averaging over all neutrino species, the total asymmetry in neutrino 
flux is of order a <~ 0.2e a b s , and the resulting kick velocity Vidck ~ 50 Bi 5 km s _1 . 
There is probably a factor of 5 uncertainty in this estimate. To firm up this 
estimate requires solving the neutrino transport equations in the presence of 
parity violation for realistic proto-neutron stars. 

4.2 Effect of Asymmetric Field Topology 

A different kick mechanism relies on the asymmetric magnetic field distribu- 
tion in proto-neutron stars (see Bisnovatyi-Kogan 1993; however, he considered 
neutron decay, which is not directly relevant for neutrino emission from proto- 
neutron stars). Since the cross section for v e (D e ) absorption on neutrons (pro- 
tons) depends on the local magnetic field strength due to the quantization of 
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energy levels for the e~ (e+) produced in the final state, the local neutrino fluxes 
emerged from different regions of the stellar surface are different. Calculations 
indicate that to generate a kick velocity of ~ 300 km s" 1 using this mechanism 
alone would require that the difference in the field strengths at the two opposite 
poles of the star be at least 10 16 G (Lai & Qian 1998b). Note that unlike the kick 
due to parity violation (see §4.1), this mechanism does not require the magnetic 
field to be ordered, i.e., only the magnitude of the field matters. 

4.3 Dynamical Effect of Magnetic Fields 

A superstrong magnetic field may also play a dynamical role in the proto-neutron 
star. For example, it has been suggested that a locally strong magnetic field can 
induce "dark spots" (where the neutrino flux is lower than average) on the stellar 
surface by suppressing neutrino-driven convection (Duncan & Thompson 1992). 
While it is difficult to quantify the kick velocity resulting from an asymmetric 
distribution of dark spots, order-of-magnitude estimate indicates that a local 
magnetic field of at least 10 15 G is needed for this effect to be of importance. 
Much work remains to be done to quantify the magnetic effects (especially when 
coupled with rotation) on the dynamics of the proto-neutron star and the super- 
nova explosion (see, e.g., LcBlanc & Wilson 1970; Thompson & Duncan 1993). 

4.4 Exotic Neutrino Physics 

There have also been several ideas of pulsar kicks which rely on nonstandard 
neutrino physics. It was suggested (Kusenko & Segre 1996) that asymmetric 
v T emission could result from the Mikheyev-Smirnov-Wolfcnstcin flavor trans- 
formation between v T and v e inside a magnetized proto-neutron star because 
a magnetic field changes the resonance condition for the flavor transformation. 
This mechanism requires neutrino mass of order 100 cV. A similar idea (Akhme- 
dov et al. 1997) relies on both the neutrino mass and the neutrino magnetic 
moment to facilitate the flavor transformation (resonant neutrino spin-flavor 
precession; see also Grasso et al. 1998; Nardi & Zuluaga 2000). More detailed 
analysis of neutrino transport (Janka & Raffelt 1998), however, indicates that 
even with favorable neutrino parameters (such as mass and magnetic moment) 
for neutrino oscillation, the induced pulsar kick is much smaller than previously 
estimated (i.e., B >• 10 15 G is required to obtain a 100 km s _1 kick). 

It is clear that all the kick mechanisms discussed in this section (§4) are of 
relevance only for B > 10 15 G. While recent observations have lent strong sup- 
port that some neutron stars ( "magnetars" ) are born with such a superstrong 
magnetic field (e.g., Thompson & Duncan 1993; Vasisht & Gotthelf 1997; Kou- 
veliotou et al. 1998,1999; Thompson 2000b), it is not clear (perhaps unlikely) 
that ordinary radio pulsars (for which large velocities have been measured) had 
initial magnetic fields of such magnitude (see als §6). 
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5 Electromagnetically Drievn Kicks 



Harrison & Tademaru (1975) show that electromagnetic (EM) radiation from 
an off-centered rotating magnetic dipole imparts a kick to the pulsar along its 
spin axis. The kick is attained on the initial spindown timescale of the pulsar 
(i.e., this really is a gradual acceleration), and comes at the expense of the spin 
kinetic energy. We (Lai, Chernoff & Cordes 2001) have reexamined this effect 
and found that the force on the pulsar due to asymmetric EM radiation is larger 
than the original Harrison & Tademaru expression by a factor of four. If the 
dipole is displaced by a distance s from the rotation axis, and has components 
/x p , yit z (in cylindrical coordinates), the force is given by (to leading order in 
f2s/c) 

15 ' 
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(2) 



(The sign is such that negative F implies Vkick parallel to the spin fl. 
dominant terms for the spindown luminosity give 
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For a "typical" situation, ji p ~ [j,^ ~ fj, z , the asymmetry parameter e = 
F/{L/c) is of order 0A((7s/c). For a given H, the maximum e max = y/OA — 0.63 
is achieved for fi p /fJ, z = and fi^/fiz — V0.4 {Qs/c). From MV = e{L/c) = 
—e(Ifi£2)/c, we obtain the kick velocity 
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where R = 10i?io km is the neutron star radius, Vi is the initial spin frequency, 
v is the current spin frequency of the pulsar, and e = (Qf — Q 2 ) ~ 1 J e dQ 2 . 
For the "optimal" condition, with fj, p = 0, fi^/fi z — VOA (/2jS / c) , and e = 
V0A [2f2 t f2/(f2 2 + Q 2 )] , we find 
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Thus, if the NS was born rotating at Vi > 1 kHz, it is possible, in principle, to 
generate spin-aligned kick of a few hundreds km s _1 or even 1000 km s _1 . 

Equations (^) and (j|) assume that the rotational energy of the pulsar en- 
tirely goes to electromagnetic radiation. Recent work has shown that a rapidly 
rotating (y > 100 Hz) NS can potentially lose significant angular momentum to 
gravitational waves generated by unstable r-mode oscillations (e.g., Andersson 
1998; Lindblom, Owen & Morsink 1998; Owen et al. 1998; Andersson, Kokko- 
tas & Schutz 1999; Ho & Lai 2000). If gravitational radiation carries away the 
rotational energy of the NS faster than the EM radiation does, then the electro- 
magnetic rocket effect will be much diminished (Gravitational radiation can also 
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carry away linear momentum, but the effect for a NS is negligible). In the linear 
regime, the r-mode amplitude a ~ £,/R (where £ is the surface Lagrangian dis- 
placement; see the references cited above for more precise definition of a) grows 
due to gravitational radiation reaction on a timescale i gr0 w — 19 (f/1 kHz) -6 s. 
Starting from an initial amplitude on -C 1, the mode grows to a saturation level 
a sa t in time t grow ln(a sat /ai) during which very little rotational energy is lost. 
After saturation, the NS spins down due to gravitational radiation on a timescale 
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(Owen et al. 1998). By contrast, the spindown time due to EM radiation alone 
is 

(7) 
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where is the surface dipole magnetic field in units of 10 13 G. Including 
gravitational radiation, the kick velocity becomes 
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where in the second equality we have replaced e by constant mean value e, and 
is defined by 
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For /3 C 1, equation (||) becomes eq. ([|); for /3 ^> 1, the kick is reduced by a 
factor 1/0. 

Clearly, for the EM rocket to be viable as a kick mechanism at all requires 
0<1. The value of a sa t is unknown. Analogy with secularly unstable bar-mode in 
a Maclaurin spheroid implies that a sa t ~ 1 is possible (e.g., Lai & Shapiro 1995). 
It has been suggested that turbulent dissipation in the boundary layer near the 
crust (if it exists early in the NS's history) may limit a sa t to a small value of 
order 10~ 2 -10~ 3 (Wu, Matzner & Arras 2000). The theoretical situation is not 
clear at this point (see Lindblom et al. 2000 for recent simulations of nonlinear 
r-modes) . 



6 Astrophysical Constraints on Kick Mechanisms 

In §§3-5 we have focused on the physics of different kick mechanisms. All these 
mechanisms still have intrinsic physics uncertainties and require more theoretical 
work. For example: 

(1) For the hydrodynamical driven kicks, one needs to better understand 
the structure of pre-SN core in order to determine whether overstable g-modes 
can grow to large amplitudes; more simulation would be useful to pin down the 
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precise relationship between the magnitudes of the initial asymmetry and the 
kick velocity; 

(2) For the neutrino-magnetic field driven kicks, more elaborate neutrino 
transport calculation is necessary to determine (to within a factor of 2) the 
value of B needed to generate (say) Vkick = 300 km s -1 ; 

(3) For the electromagnetically driven kicks, the effect of gravitational radi- 
ation (especially the r-mode amplitude) needs to be better understood. 

We now discuss some of the astrophysical/observational constraints on the 
kick mechanisms. 

6.1 Initial Magnetic Field of NS 

The neutrino-magnetic field driven kicks (§4) require initial _B>10 15 G to be of in- 
terest. While magnetars may have such superstrong magnetic fields at birth (e.g., 
Thompson & Duncan 1993; Kouvcliotou ct al. 1998,1999; Thompson 2000b), the 
situation is not clear for ordinary radio pulsars, whose currently measured mag- 
netic fields are of order 10 12 G. It is difficult for an initial large-scale 10 15 G to 
decay (via Ohmic diffusion or ambipolar diffusion) to the canonical 10 12 G on the 
relevant timescale of 10 3 — 10 7 years (see Goldreich & Reisenegger 1992). How- 
ever, one cannot rule out the possibility that in the proto-neutron star phase, a 
convection-initiated dynamo generates a "transient" superstrong magnetic field, 
lasting a few seconds, and then the field gets destroyed by "anti-dynamo" as the 
convection ceases. Obviously if we can be sure that this is not possible, then we 
can discard the mechanisms discussed in §4. 

6.2 Initial Spin of NS 

To produce sufficient velocity, the electromagnetic rocket effect (§5) requires 
the NS initial spin period Pj to be less than 1 — 2 ms. It is widely thought 
that radio pulsars are not born with such a rapid spin, but rather with a more 
modest Pi = 0.02 — 0.5 s (e.g., Lorimcr ct al. 1993). The strongest argument 
for this comes for the energetics of pulsar nebulae (particularly Crab) . But this 
is not without uncertainties. For example, a recent analysis of the energetics of 
the Crab Nebula suggests an initial spin period ~ 3 — 5 ms followed by fast 
spindown on a time scale of 30 yr (Atoyan 1999). As for the Vela pulsar, the 
energetics of the remnant do not yield an unambiguous constraint on the initial 
spin. Also, the recent discovery of the 16 ms X-ray pulsar (PSR J0537-6910) 
associated with the Crab-like supernova remnant N157B (Marshall et al. 1998) 
in the Large Magellanic Cloud implies that at least some NSs are born with spin 
periods in the millisecond range. So at this point it may be prudent to consider 
Pi ~ 1 ms as a possibility (see also §6.5). 

6.3 Natal vs. Post-Natal Kicks 

There is a qualitative difference between natal kicks (including the hydrody- 
namical driven and neutrino-magnetic field driven kicks) and post-natal kicks. 
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Because it is a slow process, the Harrison- Tademaru effect may have difficulty in 
explaining some of the characteristics of NS binaries (even if the physics issues 
discussed in §5 work out to give a large Vkick), such as the spin-orbit misalign- 
ment in PSR J0045-7319 (Kaspi et al. 1996) and PSR 1913+16 (Kramer 1999) 
needed to produce the observed precessions. For example, in the case of PSR 
J0045-7319 - B star binary: if we assume that the orbital angular momentum of 
the presupernova binary is aligned with the spin of the B star, then the current 
spin-orbit misalignment can only be explained by a fast kick with Tkick less than 
the post-explosion orbital period P rb- Similarly, a slow kick (with Tkick ^ -Forb) 
may be inconsistent with the NS binary populations (e.g., Dewey & Cordes 1987; 
Fryer & Kalogera 1997; Fryer et al. 1998). However, note that Tkick ~ t em for the 
Harrison- Tademaru effect depends on value of B [see eq. (0)], thus can be made 
much smaller than P or b (which typically ranges from hours to several months or 
a few years at most for relevant binaries) if B is large. 



6.4 Correlations Between Velocity and Other Properties of NS ? 

Despite some earlier claims to the contrary, statistical studies of pulsar popu- 
lation have revealed no correlation between Vkick and magnetic field strength, 
or correlation between the kick direction and the spin axis (e.g., Lorimer et 
al. 1995; Cordes & Chernoff 1998; Deshpande et al. 1999). Given the large sys- 
tematic uncertainties, the statistical results, by themselves, cannot reliably con- 
strain any kick mechanism. For example, the magnetic field strengths required 
for the neutrino-driven mechanisms are >10 15 G, much larger than the currently 
inferred dipolar surface fields of typical radio pulsars; there are large uncertain- 
ties in using the polarization angle to determine the pulsar spin axis; differential 
galactic rotation is important for distant NSs and cannot be accounted for unless 
the distance is known accurately and the NS has not moved far from its birth lo- 
cation; several different mechanisms (including binary breakup) may contribute 
to the observed NS velocities (see Lai, Chernoff & Cordes 2001). 

Recent observations of the Vela pulsar and the surrounding compact X-ray 
nebula with the Chandra X-ray Observatory reveal a two sided asymmetric jet 
at a position angle coinciding with the position angle of the pulsar's proper 



motion (Pavlov et al. 2000; see http: / /chandra. harvard.edu/photo/cyclcl /vela/ 



for image) The symmetric morphology of the nebula with respect to the jet 
direction strongly suggests that the jet is along the pulsar's spin axis. Analysis 
of the polarization angle of Vela's radio emission corroborates this interpretation. 
Similar evidence for spin- velocity alignment also exists for the Crab pulsar. Thus, 
while statistical analysis of pulsar population neither support nor rule out any 
spin-kick correlation, at least for the Vela and Crab pulsars, the proper motion 
and the spin axis appear to be aligned. Interestingly, both Crab and Vela pulsars 
have relatively small transverse velocities (of order 100 km s -1 ). 
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6.5 The Effect of Rotation and Spin-Kick Alignment? 

The apparent alignment between the spin axis and proper motion for the Crab 
and Vela pulsar raises an interesting question: Under what conditions is the spin- 
kick alignment expected for different kick mechanisms? Let us look at the three 
classes of mechanisms discussed in §§3-5 (see Lai, Chernoff & Cordes 2001). 

(1) Electromagnetically Driven Kicks: The spin-kick slignment is naturally 
produced. (Again, note that Pi <~ 1 — 2 ms is required to generate sufficiently 
large V kick ). 

(2) Neutrino-Magnetic Field Driven Kicks: The kick is imparted to the NS 
near the neutrinosphere (at 10's of km) on the neutrino diffusion time, Tkick ~ 
10 seconds. As long as the initial spin period Pi is much less than a few seconds, 
spin-kick alignment is naturally expected. 

(3) Hydrodynamically Driven Kicks: The low-order g-modes trapped in the 
presupernova core (M ~ 1.4M©, R ~ 1500 km) have periods of 1-2 seconds, 
much shorter than the rotation period of the core (unless the core possesses a 
dynamically important angular momentum after collapse), thus the g-modes are 
not affected by rotation. Also, since the rotational speed of the core is typically 
less than the speed of convective eddies (~ 1000-2000 km s _1 , about 20% of the 
sound speed) in the burning shell surrounding the iron core, rotation should not 
significantly affect the shell convection either. Thus the development of large- 
scale presupernova (dipolar) asymmetry is not influenced by the core rotation. 
But even though the primary thrust to the NS (upon core collapse) does not 
depend on spin, the net kick will be affected by rotational averaging if the asym- 
metry pattern (near the shock breakout) rotates with the matter at a period 
shorter than the kick timescale Tkick- Here the situation is more complicated be- 
cause the primary kick to the NS is imparted at a large radius, r s h oc k ^ 100 km 
(since this the location of the stalled shock) . To obtain effective spin averaging, 
we require the rotation period at r s h oc k to be shorter than Tkick ~ 100 ms (this 
Tkick is the shock travel time at speed of 10 4 kms^ 1 across ~ 1000 km, the radius 
of the mass cut enclosing 1.4M ). Assuming angular momentum conservation, 
this translates into the requirement that the final NS spin period Pj < 1 ms. We 
thus conclude that if rotation is dynamically unimportant for the core collapse 
and explosion (corresponding to P s » 1 ms), then rotational averaging is ineffi- 
cient and the hydrodynamical mechanism does not produce spin-kick alignment. 

The d ibove is based on the standard picture of core-collapse su- 

pernovae, which is valid as long as rotation does not play a dynamically impor- 
tant role (other than rotational averaging) in the supernova. If, on the other 
hand, rotation is dynamically important, the basic collapse and explosion may 
be qualitatively different (e.g., core bounce may occur at subnuclear density, the 
explosion is weaker and takes the form of two-sided jets; see Monchmeyer et 
al. 1991; Rampp, Muller & Ruffert 1998; Khokhlov et al. 1999; Fryer & Heger 
1999). The possibility of a kick in such systems has not been studied, but it is 
conceivable that an asymmetric dipolar perturbation may be coupled to rotation, 
thus producing spin-kick alignment. 
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It has been suggested that the presupernova core has negligible angular mo- 
mentum and the pulsar spin may be generated by off-centered kicks when the 
NS forms (Spruit & Phinney 1998). It is certainly true that even with zero prec- 
ollapse angular momentum, some rotation can be produced in the proto-neutron 
star (Burrows et al. 1995 reported a rotation period of order a second generated 
by stochastic torques in their 2D simulations of supernova explosions) , although 
Pi < 30 ms seems difficult to get. In this picture, the spin will generally be per- 
pendicular to the velocity; aligned spin-kick may be possible if the kick is the 
result of many small thrusts which are appropriately oriented (Spruit & Phin- 
ney 1998) — this might apply if small-scale convection were responsible for the 
kick. But as discussed in §3, numerical simulations indicate that such convec- 
tion alone does not produce kicks of sufficient amplitude. Therefore, spin-kick 
alignment requires that the proto-neutron star have a "primordial" rotation (i.e., 
with angular momentum coming from the presupernova core). 

Clearly, if spin-kick alignment is a generic feature for all NSs, it can provide 
strong constraints on the kick mechanisms, supernova explosion mechanisms, as 
well as initial conditions of NSs. 
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